function [ secondMembre ] = second_membre_lineique( donnees,maillage )
% second_membre_lineique retourne le second membre dû au chargement linéique pour la resolution
% numerique du calcul

a = donnees.forceLineiqueA ;
b = donnees.forceLineiqueB ;
c = donnees.forceLineiqueC ;

A = zeros(2*length(maillage),length(maillage));
Aelem = zeros(4,2) ;

for i=1:(length(maillage)-1)
   
    iCourant = (i-1)*2+1;
    jCourant = i ;
    elementCourant = maillage(i:i+1);
    
    l = elementCourant(2) - elementCourant(1) ;
    
    Aelem = [7/20 3/20 ; 1/20*l 1/30*l ; 3/20 7/20 ; -1/30*l -1/20*l]*l ;
    
    A( iCourant:iCourant+3 ,jCourant:jCourant+1 ) = A( iCourant:iCourant+3,jCourant:jCourant+1 ) + Aelem;

end

q = zeros(length(maillage),1) ;

for i=1:(length(q))
    q(i) = a*maillage(i)^2 + b*maillage(i) + c ;
end

secondMembre = A*q ;

end

